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The relaxors Pb(Zni/3Nb2/3)03 (PZN) and Pb(Mgi/3Nb2/3)03 (PMN) have very similar prop- 
erties based on the dielectric response around the critical temperature (defined by the structural 
transition under the application of an electric field). It has been widely believed that these mate- 
rials are quite different below with the unit cell of PMN remaining cubic while in PZN the low 
temperature unit cell is rhombohedral in shape. However, this has been clarified by recent high- 
energy x-ray studies which have shown that PZN is rhombohedral only in the skin while the shape 
of the unit cell in the bulk is nearly cubic. In this study we have performed both neutron elastic 
and inelastic scattering to show that the temperature dependence of both the diffuse and phonon 
scattering in PZN and PMN is very similar. Both compounds show a nearly identical recovery of 
the soft optic mode and a broadening of the acoustic mode below Tc- The diffuse scattering in PZN 
is suggestive of an onset at the high temperature Burns temperature similar to that in PMN. In 
contrast to PMN, we observe a broadening of the Bragg peaks in both the longitudinal and trans- 
verse directions below Tc. We reconcile this additional broadening, not observed in PMN, in terms 
of structural inhomogeneity in PZN. Based on the strong similarities between PMN and PZN, we 
suggest that both materials belong to the same universality class and discuss the relaxor transition 
in terms of the three-dimensional Heisenberg model with cubic anisotropy in a random field. 



PACS numbers: 

I. INTRODUCTION 

The relaxor materials of the chemical form PbBOs 
have generated much interest recently due to their large 
piezoelectric constants which are an order of magnitude 
larger than those of conventional ferroelectrics.-^'^ These 
materials are characterized by quenched disorder on the 
B site and display a diffuse phase transition with a broad 
and frequency dependent dielectric response peaked at 
the temperature T = T^ax (reproduced in the upper 
panel of Fig. Early studies of the refractive in- 

dex showed that regions of local ferroelectric order are 
formed in a paraelectric background at a temperature 
T = Td (denoted as the Burns temperature) The na- 
ture of this diffuse transition and the low temperature 
ground state has been the focus of many recent stud- 
ies. Most recently, much attention has centered on the 
material Pb(Mgi/3Nb2/3)03 (PMN) with a relatively ac- 
cessible Burns temperature of Td ~ 620 K. 

High temperature neutron inelastic scattering was 
originally investigated by Naberezhnov et alP and was 
followed by Gehring et al& who clearly showed the pres- 
ence of a soft ferroelectric mode in PMN which becomes 
overdamped below TdJ Diffuse scattering studied by Hi- 
rota et alA showed that the growth of diffuse scattering 
could be well described by the formation of phase-shifted 
polar nanoregions in which the atomic displacement was 
described by two components - a center of mass con- 
serving component and another non-conserving compo- 
nent associated with a uniform phase shift. The con- 



nection between the diffuse and phonon scattering was 
later established through a detailed phonon stud}i^ which 
showed the soft optic and the transverse acoustic modes 
to be strongly coupled. This coupling provides a nat- 
ural explanation for the presence of phase-shifted polar 
nanoregions as the center of mass conserving atomic shift 
can be attributed to the condensation of the soft optic 
mode while the non-conserving component corresponds 
to a softening of the acoustic mode. 

Unpoled PMN shows essentially no structural change 
through Tc = 213 K with only a small transverse broad- 
ening in the Bragg peaks. ^ The transition temperature 
Tc in PMN is defined by the structural transition af- 
ter poling the sample under an electric field. Despite 
no structural phase transition the soft optic mode in un- 
poled PMN was found to recover below Tc and the energy 
squared was observed to increase linearly with decreas- 
ing temperature.!?^ This is extremely surprising as such 
a recovery has previously been thought to be associated 
with a well defined structural transitiouii^ The recovery 
of the soft mode suggests the presence of a ferroelectric 
distortion in PMN, despite the absence of any change in 
the average unit cell shape. 

Until very recently, Pb(Zni/3Nb2/3)03 (PZN) and 
PMN were thought to be very different with PMN re- 
maining cubic at all temperatures and PZN undergoing 
a structural phase transition to a rhombohedral phasoi^ 
at a very well defined critical temperature Tc = 410 K. 
The Bragg scattering around Tc in both these compounds 
displays contrasting behavior as illustrated in the lower 
panel of Fig. ^ which shows a large increase in the neu- 
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FIG. 1: Schematic figure showing the temperature depen- 
dence of the dielectric constant (data taken from Ye et al.) 
and the Bragg peak intensity for both PMN (220) and PZN 
(300) which characterize rhombohedral distortion. The up- 
per figure shows the dielectric constant for 100 Hz (peaked 
at lower temperatures) and 1 MHz (peaked at higher tem- 
peratures). The dotted lines indicate the critical temperature 
measured under the application of an electric field. 



tron Bragg peak intensity for PZN at Tc but no change in 
PMN over a broad range in temperature. The structural 
transition in PZN has been reexamined in detail by a re- 
cent diffraction study using ^ 9 keV x-raysi^. Studies of 
the diffuse scattering in PZN by La-Orauttapong et al. 
showed that the diffuse scattering started near Tc, well 
below Td where the onset of diffuse scattering occurs in 
PMN, highlighting another difference between PMN and 
PZNpi^ 

The notion that the unit cell in PZN undergoes a tran- 
sition from a cubic to rhombohedral in shape at a well 
defined transition temperature Tc changed with a recent 
high-energy and penetrating x-ray study by Xu et al}^ 
and with a neutron diffraction study on PbTiOs (PT) 
doped PZN by Ohwada et alm^. Both these studies found 
that the low temperature structural phase was not rhom- 
bohedral but rather was characterized by a cubic unit 
cell, denoted as phase X. The 67 keV high energy x- 



ray study revealed phase X to be characterized by res- 
olution limited Bragg peaks suggesting a well ordered 
low-temperature structure. The connection to the PMN 
system was made recently by a high resolution neutron 
diffraction study on 10 % PT doped PMN^^ which found 
the low temperature ground state to also be character- 
ized by a cubic shaped unit cell with nearly resolution 
limited Bragg peaks. These recent studies on PZN and 
PT doped PMN show that PZN and PMN may not be 
as different as once thought. These results necessitate a 
more detailed study of the dynamical and static prop- 
erties of pure PZN and a comparison between PZN and 
PMN. 

To investigate the analogy between PMN and PZN we 
have conducted a detailed study of the structural and 
dynamical properties of PZN around the critical temper- 
ature Tc using neutron inelastic and elastic scattering. 
We will define the similarities between PMN and PZN 
through both the inelastic scattering from the phonons 
and the elastic diffuse scattering. A key difference be- 
tween these relaxor systems illustrated by the Bragg scat- 
tering will be discussed and reconciled in terms of struc- 
tural inhomogeneity in PZN. The strong similarity be- 
tween these two materials suggests they belong to the 
same universality class. We propose that the relaxor 
transitions can be understood in terms of the three- 
dimensional Heisenberg model with cubic anisotropy in 
the presence of a random field. This model is very attrac- 
tive as it unifies the many temperature scales measured 
by various experimental probes. 



II. EXPERIMENTAL DETAILS 

Neutron scattering experiments were performed at the 
NRU reactor, Chalk River Laboratories. Elastic scatter- 
ing studies of the diffuse scattering and the Bragg peaks 
were conducted on the N5 spectrometer by holding the 
incident and final energies fixed at 14.8 meV. Diffuse 
scattering results were obtained using a flat Graphite 
(002) monochromator and analyzer. Horizontal coUima- 
tion was set by SoUer slits from the reactor to detector 
at 77'-33'-5-35'-60' (where S denotes the sample) and a 
Graphite filter was used on the scattered side to remove 
higher order neutrons. For elastic studies of the Bragg 
peaks we have matched the d-spacings of the sample to 
that of the (220) reflection of a Germanium analyzer with 
a small mosaic of ^3'. By matching the d-spacings of the 
sample and analyzer crystal we have obtained good res- 
olution allowing us to study the small broadening of the 
Bragg peaks below Tc. For elastic measurements of the 
Bragg peaks we have also set the horizontal collimation 
to 77'-6'-5-12'-60'. 

Inelastic measurements studying both the transverse 
acoustic and optic modes were conducted on the C5 spec- 
trometer by fixing the final energy to 14.5 meV and vary- 
ing the incident energy. A variable focusing graphite 
(002) monochromator and a flat graphite analyzer were 



used with the horizontal coUimation set at 33'-29'-5'-48'- 
72' to study the acoustic mode and 33'-48'-5'-51'-72' to 
study the optic mode. A graphite filter was also used on 
the scattered side to filter out higher order neutrons and 
a cold sapphire filter was used before the monochromator 
to reduce the fast neutron background. 

Single crystals of PZN were grown by spontaneous nu- 
cleation from high temperature solution using PbO as 
flux, based on the technique previously describedjiSc A 
single crystal with a volume of 0.2 cc and a weight of 1.8 
g was selected for the experiments. It was oriented with 
the largest face parallel to the (100) cubic plane. The 
sample of PZN was mounted, using tantalum wire, to a 
boron-nitride post with a copper insert to ensure good 
thermal contact. The sample was mounted in an orange 
cryofurnace such that reflections of the form (HKO) lay 
in the scattering plane. The room temperature lattice 
constant was measured to be 4.04 A and therefore 1 re- 
ciprocal lattice unit (r.l.u.) corresponds to 27r/a = 1.55 
A~^. Because PZN decomposes at very high tempera- 
tures we have kept our high temperature measurements 
below 550 K to avoid any possible damage. Unfortu- 
nately, this means that the studies above or around the 
Burns temperature Td, which exceeds 670 K, are not pos- 
sible in PZN. 



III. BRAGG SCATTERING 

To investigate the Bragg peak profile as a function of 
temperature around Tc we have used a Germanium an- 
alyzer with a small mosaic spread of ^ 3'. Based on 
low-energy x-ray results showing a transition from cubic 
to rhombohedral we conducted scans around the (110) 
position. If a rhombohedral phase exists at low tem- 
peratures the Bragg peak should split along the longitu- 
dinal [110] direction only. Instead of such a distortion 
we observe a broadening of the (110) Bragg peak along 
both the longitudinal [110] and transverse [110] direc- 
tions (Fig. O. This is in contrast to the behaviour in 
PMN which shows only transverse broadening below Tc 
(Ref. 9). We have also investigated the temperature de- 
pendence of the (200) Bragg peak which is found to also 
display both transverse and longitudinal broadening be- 
low Tc (Fig. ISJ. As mentioned in the next section, we 
find no observable temperature dependent diffuse scat- 
tering around the (200) position. This means that the 
broadening of the Bragg peaks shown in Fig.Ocannot be 
associated with the diffuse component. 

These results contrast with the results of both low- 
energy and high energy x-ray measurements. The ~ 9 
keV x-ray study by Lebon et al. showed the ground state 
of PZN at low temperatures to have a rhombohedral unit 
cell. The high-energy 67 keV x-ray results by Xu et al. 
found the unit cell to be nearly cubic and to be very well 
ordered characterized by resolution limited Bragg peaks. 
Xu et al. also confirmed the result of Lebon et al. by 
showing that a splitting of the (111) Bragg peak charac- 
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FIG. 2: Contour plots of the (110) reflection as scanned in the 
(HKO) plane above (a) and below (6) the critical temperature 
Tc. These data were taken using a Germanium analyzer with 
a small mosaic spread. Dashed lines indicate longitudinal 
[110] and transverse [110] directions. 

teristic of a rhombohedral distortion existed when a lower 
incident beam energy of 32 keV was used. The contrast- 
ing results presented by x-ray scattering with various in- 
cident energies can be reconciled based on the relative 
penetration depths. As noted by Xu et at, the lower en- 
ergy x-rays in the range 10-32 keV have a penetration 
depth in PZN of only ~ 10-60 fim whereas x-rays with 
an incident energy of 67 keV have a penetration depth 
of ^ 400 fim which should give a more accurate picture 
of the bulk properties. Therefore, in PZN rhombohedral 
order is formed in the skin of the sample whereas the 
bulk has a well ordered pseudocubic structure. This new 
bulk ground state has been referred to as Phase-X. 

Our neutron results showing both transverse and lon- 
gitudinal broadening of the Bragg peaks are consistent 
with the picture obtained from x-rays. Since neutrons 
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FIG. 3: Plots of the transverse, longitudinal widths, and 
the integrated intensities as a function of temperature. The 
widths were obtained from Lorentzian fits to the observed 
scattering. In all cases there is clear anomaly at the critical 
temperature Tc. The solid lines are guides to the eyes. 



probe the entire sample volume the observed Bragg peak 
profile is a combination of both the bulk and the near- 
surface region. The presence of both transverse and lon- 
gitudinal broadening is consistent with finite size effects 
resulting from domains. Thus, based on the x-ray re- 
sults, the broadening along both directions is very likely 
the result of structural inhomogeneity in the bulk sam- 
ple. The size of these domains can be estimated from the 
half- width (estimated from a Lorentzian profile) of the 
Bragg peaks at low temperatures to be ~ 200-400 A. This 
estimate is consistent with the value of 700 A estimated 
by Xu et al. based on the low energy x-ray results. 

The existence of finite size effects resulting from struc- 
tural inhomogeneity is also consistent with the tempera- 
ture dependence of the integrated intensity displayed in 
Fig. 13 for both the (110) and (200) Bragg peaks. Both 
Bragg peaks show a large increase in the integrated in- 
tensity at Tc characteristic of a release of extinction at 
low temperatures. Such an effect is the result of a broad- 
ening of the overall mosaic spread of the crystal allowing 
better penetration of the neutron beam into and out of 
the crystal. Detailed analysis of extinction and its de- 
pendence on the domain size is discussed elsewhereiSi 

Thus, our results are consistent with the high-energy 
x-ray results showing the absence of rhombohedral order 
at low temperatures. We do observe a longitudinal line 
broadening indicating structural inhomogeneity. In the 
sense that no rhombohedral order is formed in the bulk 
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FIG. 4: Diffuse scattering around the (300) position at several 
temperatures. The diffuse scattering is clearly present well 
above the critical temperature Tc = 410 K and suggestive 
that diffuse scattering starts at the Burns temperature Td- 
The inset shows the unconvoluted HWHM of the (300) Bragg 
peak as a function of temperature. 



at low temperatures, the structural properties of PZN 
are similar to those of PMN. 



IV. DIFFUSE SCATTERING 

The absence of a clear rhombohedral distortion at 
low temperatures necessitates a revaluation of the dif- 
fuse scattering. Previous diffuse scattering results by 
D. La-Orauttapong et al^ have indicated that the tem- 
perature dependent diffuse scattering has an onset near 
Tc ^ 410 K. This result contrasts with that in PMN 
where the diffuse scattering starts at the much higher 
temperature of Td where polar nanoregions are formed. 

We have studied the diffuse scattering and its tem- 
perature dependence around the (300) reffection where 
the diffuse scattering is strong. We have found that the 
lineshape of the diffuse scattering at low temperatures is 
qualitatively similar to that of PMN studied in detail by 
You et al?"^ These authors have found that the diffuse 
intensity in PMN is elongated along both the [110] and 
[110] directions at the (300) position. Figure 01 shows the 
temperature variation of the (300) diffuse profile indicat- 
ing the presence of diffuse scattering even above 500 K. 
The presence of a diffuse component above Tc is clearly 
demonstrated in Fig. [S] which shows the temperature de- 
pendence of the diffuse intensity at different g-positions 
with a constant background indicated by a dashed line. 
It is again clear that diffuse scattering is still present 
at 550 K. These results indicate that diffuse scattering 
starts well above Tc and is suggestive of an onset at the 
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FIG. 5: Measured intensity at (3, q, 0) as a function of 
temperature for various values of q. Large values of q show a 
gradual and continuous growth with decreasing temperature 
of the diffuse scattering through Tc while smaller values of 
q show evidence of critical scattering. The horizontal dotted 
line indicates the constant background. 



Burns temperature Td similar to the case in PMN. As 
also displayed in Fig. 0] this continuous growth of the 
diffuse scattering is in contrast to the transverse HWHM 
of the Bragg peak which shows a sharp anomaly at Tc. 
The gradual growth of the diffuse intensity starting at 
temperatures in excess of 500 K differs with the sugges- 
tion that the diffuse scattering starts at ^ 450 and 
is consistent with previous diffuse work on PMN which 
clearly showed the onset of a diffuse component at the 
Burns temperature.^'* 

Another important characteristic of diffuse scattering 
in PZN is shown in Fig.|Sl For large values of q the diffuse 
scattering grows continuously through Tc but at smaller 
values of g a clear peak is observed in the measured inten- 
sity suggesting the presence of critical scattering. Even 
though this is consistent with the diffuse scattering in 
PMN-PT^^ these results differ from the current under- 
standing of diffuse scattering in pure PMN which shows 
a continuous change through Tc.'^'^ Therefore the diffuse 
scattering around Tc could be interpreted as consisting 
of scattering from the polar nano-regions formed at Td 
and critical scattering from some component of the sam- 
ple undergoing a structural phase transition at Tc (see 
Fig0. Based on our analysis of the Bragg scattering 
this extra component of the sample which undergoes a 
structural transition is located in the near-surface region 
and is the same component studied by Lebon et at using 
low energy x-rays. 

We believe that the structural inhomogeneity is the un- 
derlying reason for the apparent discrepancy between our 



diffuse scattering results and those of a previous diffuse 
study on pure unpoled PZN by D. La-Orauttapongii^ In 
that study, a two phase critical scattering model was ap- 
plied to the diffuse intensity around the (110) position 
which used a Lorentzian function to describe the diffuse 
scattering and a Gaussian to describe the Bragg peak. 
This analysis showed the onset of the diffuse Lorentzian 
component near Tc and no change in the Bragg peak 
width. We believe that this model is oversimplified as 
the scattering should be the sum of the diffuse compo- 
nent, the Bragg peak, and the critical scattering coming 
from the skin of the sample which does enter a rhombohe- 
dral phase at low temperatures. Therefore we speculate 
that the diffuse Lorentzian component in the fits by D. 
La-Orauttapong et al. reflected mostly the critical com- 
ponent not characteristic of the bulk phase (phase-X) of 
the sample. 

To test further the similarity with PMN we have 
checked for diffuse scattering at the (200) and (110) po- 
sitions at temperatures above and below Tc- We have 
found no observable change in the scattered intensity 
around the (200) Bragg peak position and only a small 
change around the (110) Bragg peak. This implies the 
absence of diffuse scattering around (200) and only a 
weak diffuse component around (110). These results for 
the structure factors are predicted by the phase shifted 
polar nanoregion model used to describe the diffuse scat- 
tering in PMN. In this model the atomic displacements 
are divided into a component which conserves the center 
of mass and another, associated with an overall phase 
shift, which does not. The component associated with 
the center of mass conserving optic modes is determined 
by a linear combination of the Last and Slater modes J4 
Based on a calculation similar to that conducted by Hi- 
rota et al. and assuming the same parameters for the 
product of the atomic displacement and the atomic mass 
at each lattice site as in PMN we find the calculated 
phase shift for PZN to be Sghift ~ 0.56 (normalized to 
the total shift of the Pb atom). This is quite similar to 
that of PMN, Sshift ~ 0.58. Hirota et al. found by choos- 
ing the correct weight of both Last and Slater modes and 
using the calculated phase shift, the intensity of the dif- 
fuse scattering at the (200) position could be made small 
compared to the diffuse scattering at the (110) and (300) 
positions. This is the same trend observed in PZN. 

The absence of any observable temperature dependent 
diffuse component around the (200) position is also im- 
portant as it illustrates that the line broadening observed 
in the Bragg scattering is coming from structural effects 
with an onset at Tc and not from a diffuse component (see 
Fig. 0). This result contrasts with the previous sugges- 
tion that the apparent change in the Bragg peak widths 
were a result of a strong diffuse temperature dependence 
around Tc.^^ Clearly, the Bragg peaks show an anomalous 
linewidth broadening beginning at the critical tempera- 
ture Tc in PZN (see Fig.EJ). 

Our analysis based on the temperature dependence 
and structure factors link the diffuse scattering directly 
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with the formation of polar nanoregions. Despite the 
fact that we are not able to reach the Burns temperature 
Td ^ 700 K our results suggest that the onset of the dif- 
fuse scattering starts at (where polar nanoregions are 
formed). Our results are strongly analogous to those in 
PMN where neutron scattering has shown the presence of 
a diffuse component at the accessible Burns temperature 
Td ^ 620 K. 



V. LATTICE DYNAMICS 

In the previous sections, we have shown that PZN has 
similar diffuse scattering characteristics to those of PMN. 
Another important feature closely related to the relaxor 
mechanism are the properties of the transverse optic and 
acoustic phonon modes. These phonon modes in PMN 
have been well studied. Striking features in PMN appear 
around Tc, such as recovery of the soft optic mode below 
Tc and disappearance of the acoustic phonon line broad- 
ening. In this section, we report the transverse optic and 
acoustic phonons in PZN measured around Tc and com- 
pare in detail these results to those of PMN to illustrate 
the similarity between these two relaxor systems. Firstly, 
we introduce the lineshape used to analyze our phonon 
results. 



A. Lineshape: Simple Harmonic Oscillator 
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FIG. 6: Scans of the transverse optic mode at (220). The solid 
line is the sum of the convoluted Lorentzian and a Gaussian 
centered at cu—Q. The dashed line is the Lorentzian contribu- 
tion and the dotted-dashed line is the Gaussian at zero energy 
transfer. 



Measurement of the neutron scattering intensity pro- 
vides a direct measure of S(Q, w), which is related to 
the imaginary part of the susceptibility x"(Q, uj) by the 
fluctuation dissipation theoremS^, 

S{q,u;)^TT-'[n{u;) + l]x"{cL,u;), (1) 

where the energy transferred to the sample is defined by 
hio = E^- Ef and n(w) = 1/(6?"^/'=^^ - 1) is the Bose 
factor. In order to obtain detailed information on the 
linewidth from the inelastic spectrum, a model must be 
convolved with the resolution function and then fit to 
the observed scattering. To analyze the phonon scatter- 
ing as a function of temperature we have taken x" to 
be described by the formula for the damped simple har- 
monic oscillator. We have used the equation given by the 
antisymmetrized linear combination of two Lorentzians, 



we have approximated the dispersion to be linear, 
'^o(q) = c|q|, where c is the phonon velocity. For op- 
tic modes, we have set the dispersion to have the form 
^o(q)^ = + (A|q|)^, where fio is the soft-mode energy 
at g = and A is a constant for a given q directioni2LS& 
The antisymmetrized Lorentzian lineshape obeys de- 
tailed balance2& which requires x" to be an odd func- 
tion in energy. By making several substitutions22, the 
lineshape in our analysis can be shown to be equivalent 
to those employed in previous studies which used a fre- 
quency independent damping constant, Fq.^'' The mea- 
sured intensity was fitted to Eq. [21 convolved with the 
instrumental resolution and a non-convoluted Gaussian 
function was used to describe the component at w=0. A 
constant was used to describe the overall background. 



B. Soft Optic Mode 



A 

where F(cj) is the frequency dependent half-width-at- 
half- maximum (HWHM), Ci^o(q) is the undamped phonon 
frequency, and A is the amplitude. For acoustic modes. 



The soft transverse optic mode was studied by scans 
with fixed Q at the (2, 2, 0) position. Figure displays 
scans over a broad range of temperature below Tc- The 
solid lines are the results of fits to the antisymmetrized 
Lorentzian lineshape convoluted with the instrumental 
resolution. The data shows a clear softening and damp- 
ening of this mode with increasing temperature until the 
optic peak is essentially unresolvable around Tc in Fig.|^ 
(d). We have also confirmed that the soft optic mode is 
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PZN : Q = (2, 2, 0), Ef = 14.5 meV, 33-48-51-72 
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FIG. 7: Transverse optic fitted parameters as a function 
of temperature. The upper panel is the phonon frequency 
Wo plotted as a function of temperature. The HWHM F is 
plotted in the lower panel. 



unobservable above Tc up to at least 550 K. 

The phonon frequency squared, {hujo)'^, and the half- 
width are plotted in Fig.^l The square of the phonon fre- 
quency, Wq, is related to the dielectric constant e via the 
Lyddane-Sachs- Teller (LST) relation which states that 
(1/e) oc (hcuo)'^- The linear relationship between the 
square of the phonon frequency and temperature has 
been established both above and below Tc in conven- 
tional ferroelectrics^ and is explained in the theory of 
Cochrani^. A linear recovery of the soft optic mode 
has been observed in PMN at low temperatures.^^ These 
combined results for both PMN and PZN suggests the ex- 
istence of a well-defined ferroelectric distortion over the 
length scale characterized by the resolution (~ 100 A). 
Therefore these results point to at least a local ferroelec- 
tric polarization in the unpoled systems. Local or mi- 
crodomain polarization in PZN has also been observed 
by Raman and optic studies'^"'^ and it was suggested that 
the unpoled system had only orientational order and not 
translational order. 

As Tc is approached from low temperatures, the op- 
tic mode becomes gradually more damped and above the 
critical temperature it becomes completely overdamped 
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FIG. 8: The soft mode position (huio)^ as a function for 
reduced temperature t = T/T^ for both PMN and PZN. The 
data for both materials overlaps extremely well illustrating 
the strong similarity between these two materials. The dashed 
line above t=l is a schematic representation of 1/e taken from 
Viehland et al. 



and is unobservable up to 550 K, the highest temper- 
ature studied. This is the "waterfall" region where the 
optic mode becomes overdamped as a result of the forma- 
tion of polar nanoregions. Polar nanoregions were first 
suggested by Burns and Dacoli based on changes in the 
index of refraction below a high temperature denoted as 
the Burns temperature T^. They found the data to be 
well described by local regions of ferroelectric order in 
a paraelectric background. In this temperature range, 
the uncorrelated polar nanoregions with random polar- 
ization prevent the propagation of the optic mode. This 
effect has been clearly observed in PMN where Td ~ 620 
K is easily accessible^ illustrating that the temperature 
dependence of the soft optic mode in PZN is qualita- 
tively similar to that in PMN. The waterfall effect and 
polar nanoregions are reviewed elsewhere in the context 
of coupled mode theory'' and has been verified by other 
groups''^ in this system. 

The recovery of the soft optic mode below Tc is found 
to be quantitatively analogous to the recovery observed 
in PMN. In Fig. |S1 we plot [fiujo)"^ as a function of the 
reduced temperature t = T/Tc for both PMN and PZN 
showing that the slope of the phonon recovery is the ex- 
actly same in both PMN and PZN. We also reproduce 
the high-temperature data for the dielectric constant in 
PMN. The data for PMN are taken from Wakimoto et 
ali^ and the data for the high-temperature dielectric con- 
stant (1/e) is taken from Viehland et alm^ Since the soft 
optic mode energy characterizes the degree of the ferro- 
electric distortion, Fig. |H1 demonstrates that both PZN 
and PMN develop an analogous ferroelectric distortion 
at low temperatures. 
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FIG. 9: Scans of the transverse acoustic (TA) mode above 
and below Tc showing a small broadening above the critical 
temperature ~ 400 K. The solid and dashed lines are the 
results of least-squares fits to the lineshape discussed in the 
text convoluted with the instrument resolution. The solid bar 
indicates the calculated instrumental resolution FWHM. 



C. Acoustic Mode 

Studies of the transverse acoustic mode (TA) were 
made by conducting scans with fixed Q at the (2, 0.12, 
0) position. Representative scans above and below Tc are 
displayed in Fig. I^where the solid and dashed curves are 
the results of least-squares fit to the lineshape previously 
described convoluted with the resolution. A comparison 
of the results above and below Tc shows a slight broad- 
ening at high temperatures. 

The half-width is plotted as a function of temperature 
in Fig. ^1 (a) and shows a sharp anomaly at Tc while 
the integrated intensity is essentially constant over the 
temperature range studied. The recovery of the trans- 
verse acoustic mode for temperatures below Tc is similar 
to that observed in PMN" and FT doped FMN^^. In 
both these cases the dampening of the acoustic mode 
appears below Td and is attributed to the formation of 
polar nanoregions. The recovery of the acoustic mode 
is associated with the development of local correlations 
between polar nanoregions resulting in local ferroelectric 
polarization. Since the acoustic modes are associated 
with center of mass motion they are less sensitive to the 
random field associated with the polar nanoregions and 
therefore should only be slightly broadened compared to 
the complete overdamping of the optic mode already pre- 
sented. 

Also, the recovery of both the optic and acoustic modes 
at Tc suggests the coupling of these modes, and therefore 
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FIG. 10: Fitted parameters as a function of temperature. 
Figure (a) shows the HWHM F as a function of temperature. 
The integrated intensity is plotted in figure (&) and has an 
essentially fiat temperature dependence. The solid lines are 
guides to the eye. 



the dynamics should be well described by coupled-mode 
theory.^* This was observed to be the case in PMN where 
the temperature and the zone dependent lattice dynam- 
ics could be well described by coupled mode theory.^ In 
FMN the appearance of diffuse scattering at Td was in- 
terpreted in terms of the condensation of the coupled soft 
mode. Even though we are not able to probe tempera- 
tures near the Burns temperature Td the strong analogy 
with PMN suggests that the broadening observed in the 
acoustic mode starts at Td and is associated with the 
condensation of the coupled soft mode, showing another 
similarity between FZN and PMN. 



VI. DISCUSSION AND CONCLUSIONS 

We have shown that PMN and PZN are essentially 
identical in terms of the temperature dependence of the 
diffuse scattering and the phonon scattering. The tem- 
perature dependence of the diffuse scattering is sugges- 
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tive of an onset not at the critical temperature Tc, but 
at the much higher Burns temperature where polar 
nanoregions are formed. The temperature dependence 
of the phonon scattering has also been shown to be ex- 
actly analogous to that in PMN. The recovery of the 
soft optic mode in PZN is identical to the observed re- 
covery in PMN. At temperatures above Tc, the optic 
mode becomes completely overdamped due to the polar 
nanoregions and the acoustic mode also becomes slightly 
damped. Both of these results are identical to those 
found in the PMN system. 

Despite these similarities, the Bragg peaks in PZN 
show the onset of both transverse and longitudinal broad- 
ening at Tc parallelled by a large growth in the integrated 
intensity. We have argued based on a direct compari- 
son of our neutron results to the previous high energy 
x-ray results of Xu et al. that this indicates the pres- 
ence of structural inhomogeneity, probably in the bulk of 
the sample. As illustrated in Fig. ^ PMN does not show 
any anomaly in the integrated intensity. No change is 
also observed in the longitudinal width and only a sub- 
tle broadening is observed in the transverse width,^ This 
suggests that the key difference between PMN and PZN 
is the presence of macroscopic structural disorder at low 
temperatures in PZN. 

A similar discrepancy in the structural properties as 
a function of incident x-ray energy has been observed in 
the doped PZN-PT system where the physical proper- 
ties have been observed to vary as a function of sample 
thickness. This has been directly shown by Noheda et al. 
in the PZN-PT system and is illustrated by several ex- 
periments using both high and low energy x-rays^^i^^*^ 
Our result, combined with those of previous x-ray studies 
conducted at various energies, provides further support 
to the growing evidence that structural inhomogeneity 
is a general feature of PZN and possibly to all relaxor 
systems. This point is further illustrated by the neu- 
tron powder diffraction work by Iwase et al. suggest- 
ing the existence of two structural phases below T^Mi^ 
We believe that previous probes sensitive to near-surface 
effects, like low energy x-rays with a small penetration 
depth and powder diffraction, showing rhombohedral or- 
der have been interpreted as a bulk result and therefore 
marking a clear difference between the PMN and PZN 
systems, whereas in fact the bulk structural phases are 
quite similar. 

The recovery of the soft optic mode at low tempera- 
tures in both PMN and PZN points to the presence of 
ferroelectric polarization on a length scale of at least ~ 
100 A (defined by the resolution) in both systems de- 
spite the clear absence of any rhombohedral order in the 
bulk phase (phase X). This suggests that the ferroelec- 
tric order parameter (the atomic shift or the polariza- 
tion) is decoupled from the unit cell shape. This same 
effect was recently observed in strained SrTiOa filmsi^ 
In this system strong superlattice peaks characteristic of 
a structural transition were observed without any change 
in the c-axis lattice parameter. This was interpreted as 



the internal degrees of freedom (in that case the TiOg 
rotations) becoming uncoupled from the overall lattice 
shape. In PZN and PMN a ferroelectric polarization is 
indicated by the recovery of the soft optic mode with- 
out the presence of a rhombohedral distortion. This fact 
suggests that the atomic shift is decoupled from the unit 
cell shape in phase X. This picture is similar to that of 
Lebon et al.^^ who have suggested that local orientational 
order (microdomains) characterizes the low temperature 
ground state (phase-X). We note though that phase-X 
is structurally well ordered as characterized by the reso- 
lution limited Bragg peaks observed by Xu et al. using 
high energy x-rays. 

Therefore, the characteristic feature of relaxors PZN 
and PMN is possibly a weak coupling between the ferro- 
electric order parameter and the structure as seen in the 
unit cell shape. In the case of conventional ferroelectric 
systems like BaTiOs the coupling is strong resulting in a 
well defined structural transition with long-range trans- 
lational and orientational order. In the case of PMN the 
coupling is weaker resulting in the development of local 
ferroelectric polarization without any change in the unit 
cell shape. In PZN, the coupling is still weaker; thus the 
bulk phase shows no change in unit cell. However the 
coupling is slightly larger than that of PMN, therefore 
resulting in a well defined change in the unit cell shape 
with the aid of some extra strain as found near the sur- 
face of a sample. 

The strong similarity between PMN and PZN and the 
nature of the relaxor phase transition in both materi- 
als, as discussed above, point to a common model for 
the relaxor transition. Hereafter we discuss a universal- 
ity class which possibly describes the relaxor transition 
from the viewpoint of an analogy of the polarization vec- 
tor to magnetic spins. When trying to formulate a model 
for the relaxor transition it is important to consider that 
there are two key temperatures which need to be de- 
scribed; the Burns temperature and the lower crit- 
ical temperature T^. Also, the model needs to explain 
the reason for the history dependence of measured quan- 
tities around T^. To satisfy these constraints, we pro- 
pose the three-dimensional Heisenberg model with cubic 
anisotropy universality class in the presence of random 
fields. Variants of this model have been considered pre- 
viously although our particular formulation and analysis 
appear to be unique42i^i*^ The mapping is as follows: 
the Heisenberg spin corresponds to the local ferroelec- 
tric polarization, the cubic anisotropy reflects the pref- 
erential orientation of the polarization along the (111) 
axes and the isotropic random magnetic field corresponds 
to the randomly oriented local electric fields originat- 
ing from the differing charges of the Nb^+ and Zn^+ 
(Mg^+). Following a suggestion by AharonyH, we con- 
sider the case where the Heisenberg term Hueisenberg 
dominates over the random field H^ip which is larger 
than the contribution from the cubic anisotropy Hcubic 

[H Heisenberg > H RF > Hcubic)^ 

In this scenario there would be two distinct tempera- 
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ture regions. At high temperatures the cubic anisotropy 
is irrelevant and therefore the system should behave like 
a Heisenberg model in a random field. In such a case the 
excitation spectrum is characterized by Goldstone modes 
and therefore no long-range order is expected in the pres- 
ence of random fields4^ The second temperature scale 
appears at low temperatures where the cubic anisotropy 
becomes relevant and therefore the system should be sim- 
ilar to that of an Ising system in the presence of a random 
field. In such a system we would expect long-range order 
to occur at equilibrium or for cooling in the presence of 
random fields, which typically yields a non-equilibrium 
state, local order with history dependent effects. 

In the actual relaxor materials PMN and PZN, the 
Burns temperature Td marks the entry into the higher 
temperature region discussed above. Therefore, since the 
order parameter has a continuous symmetry**^ and the 
excitation spectrum is characterized by gapless modes, 
the system does not order but forms polar nanoregions in 
a paraelectric background. On the other hand, the crit- 
ical temperature characterizes the lower temperature 
region where the the system behaves more like an Ising 
model in the presence of a random field. This explains 
the local ferroelectric distortion characterized by the re- 
covery of the soft-optic mode. Even though, in equilib- 
rium, a 3D Ising system should have long-range order in 
the presence of small random fields, as found in magnetic 
systems, non-equilibrium effects with long time scales be- 
come dominant. The presence of non-equilibrium effects 
may explain the lack of true long-range order at low tem- 
peratures and the history dependence of physical prop- 
erties like the linear birefringence. We also note that the 
presence of the phase-shifted polar nanoregions may also 
pose another energy barrier for the ordered phase at the 
critical temperature. 

The low-temperature properties (below Tc) of PZN 
and PMN are exactly analogous to that of Mno.5Zno.5F2 
in a magnetic field^SiH which is well described by a three- 



dimensional random field Ising model. In this system 
the bulk low temperature phase is characterized by non- 
equilibrium effects and therefore a disordered metastable 
state was observed. Importantly, an ordered phase was 
observed in the near-surface region of the sample, ex- 
actly analogous to the rhombohedral order observed in 
the near-surface region of PZN. The strong similarity be- 
tween these two systems further affirms the model we 
have proposed. Our low-temperature phonon results are 
also qualitatively consistent with the observation of well- 
defined spin-waves in Mno.5Zno.5F2 in the presence of 
random fieldsi^ A detailed, quantitative analysis of this 
model would be extremely valuable. 

We have shown that PMN and PZN are very similar 
relaxor systems through identical behavior in the diffuse 
and phonon scattering. We do observe structural inho- 
mogeneity at low temperatures in the PZN system in 
contrast to PMN. We speculate this difference is the re- 
sult of stronger coupling in PZN between the ferroelec- 
tric order parameter and the unit cell shape. We sug- 
gest a consistent explanation for the relaxor transition 
in terms of the three-dimensional Heisenberg model with 
cubic anisotropy in the presence of random fields. 
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